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Abstract 

We derive a reduction formula that expresses the survival rate for the three-flavor neutrino 
oscillation using the two-flavor one to next-to-leading order when there is one resonance due to the 
matter effect. We numerically find that the next-to-leading reduction formula is extremely accurate 
and the improvement is relevant for the precision test of solar neutrino oscillation and the indirect 
measurment of CP violation in the leptonic sector. We also derive a reduction formula, which is 
slightly different from that previously obtained, in the case when there are two resonances. We 
numerically verify that this reduction formula is quite accurate and is valid for a wider parameter 
region than are those previously obtained. 
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I. INTRODUCTION 



The physics of neutrino oscillation is currently under very active investigation, since it 
leads to physics beyond the standard model. Two-flavor neutrino oscillation, however, is 
adopted in most analyses of the data, although everyone knows that there are three active 
neutrino flavors. Two-flavor oscillation is easy to investigate in comparison with three-flavor 
oscillation because there are only two parameters: a mass-squared difference and a mixing 

ity are known in the two- 
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angle. In addition, some exact solutions of the oscillation probabi 
flavor oscillation scheme even in the presence of a matter effect 
make the analysis realistic, however, we need to work in the full three-flavor system, where 
we have six parameters, two mass-squared differences, three mixing angles, and one CP 
phase. The analysis using three-flavor oscillation is particularly important when we compare 
the outputs from different experiments, sensitive to different mass-squared differences and 
mixing angles. The importance of the analysis in the full three-flavor context for terrestrial 
neutrino oscillation experiments has been discussed jil llf|. 

The simplest way to investigate the three-flavor oscillation is by relying on a numerical 
calculation. There is no serious technical difference between the two-flavor oscillation and 
the three-flavor one in numerical calculations 35]. The parameter space, however, becomes 
very large, six dimensional, and is difficult to exhaust. What is worse, we cannot easily 
understand the physical consequences from the numerical results intuitively, even if the 
parameter space is exhausted. 

A more elegant way to investigate the three-flavor oscillation is to reduce the three-flavor 
oscillation to an effective two-flavor one using the hierarchy between the two mass differences 
or the smallness of some mixing angle. For example, the following reduction formula for the 
survival probability of the electron neutrino, S% v is found in the case that one mass difference 
is much larger than the other one and the matter effect, A 12 , A(x) -C A 13 

S 3 u = S 2v ,e f f cos 4 #13 + sin 4 #13 . ( 1 ) 

Here, #13 is the mixing angle defined in Sec. Ill Al while S^.e// is the survival rate calculated 
in the effective two-flavor scheme with effective matter effect A e ff(x) = cos 2 6i 3 A(x). We 
show the definitions of each qua ntity in detail in Sec. [ ll This relation is often used in the 



analysis of the solar neutrino 
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19| . Although it is exact in the limit 



Ai 2 /Ai 3 , A(x)/Ai 3 — > 0, the actual hierarchy is not so good. Taking the large mixing angle 



Mikheyev-Smirnov-Wolfenstein (LMA-MSW) solution for the solar neutrino problem and 
the allowed range of A13 from the analysis of the atmospheric neutrino oscillation, we find 
that the hierarchies are rather mild, A 12 /A 13 < 1/15, A/A 13 < 1/10 [36]. This indicates 
that the above relation possibly has an error around several percents or more. 

Such an error is not so important in just verifying the existence of the neutrino masses. 
Current interest in neutrino physics is, however, not only in verifying the finite masses but 
also in determining the precise values of the parameters. The allowed region of the mass 
differences and the mixing angles may be affected due to error in the formula used in the 
analysis. This error would be more serious in attempts to observe the CP violation in 
the leptonic sector usarg the sizes oX^arity triple, srnce it recces a more precise 



determination of the survival rate 
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One of the aims of the present paper is to examine how precise the reduction formula is. 
We then propose a next-to-leading order reduction formula, which is surprisingly precise. 
We assume only the mass hierarchy A12/A13 ^ 1/15 and do not impose any restriction on 
#13 in order to keep our analysis as general as possible. While we average the survival rate 
with respect to the final time, we do not average it with respect to the initial time. The 
reduction formula obtained therefore can be used without averaging over the initial time. 

Another example of the reduction formula is the following relation for two successive 

Q 



resonances 



S3* = P lPh cos 2 9 12 cos 2 13 + P H {1 - P L ) sin 2 6 X2 cos 2 6 13 + (1 - P H ) sin 2 9 13 . (2) 

Here, #12 is the mixing angle defined in Sec. Ill Al and and Ph are the jump probabilities 
for resonances with lower and higher number densities, respectively. This relation is used 
for the investigation into neutrino oscillation inside a supernova, which has two successive 
_s d ue to the high matter density MSf 1 



|22| • It is also used for investigation 



into hypothetical O(GeV) neutrinos from the sun produced by the annihilation of weakly 
interacting massive particles (WIMPs), which have two successive resonances due to the 
higher neutrino energy j30j. In Ref. 30j], an improved reduction formula Eq. (|55|) is derived. 
The reduction formulas Eq. (j2j) and Eq. (J55|) . are expressed by the jump probability, while 
the reduction formula Eq. (JJJ is expressed by the survival rate. We express the reduction 
formula by the survival rate also in this case and numerically study the difference between 
these expressions. 
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The present paper is organized as follows. In Sec. Ill Al we derive the next-to-leading order 
reduction formula in the case of one resonance with the matter effect. The main result in 
this section is the next-to-leading reduction formula Eq. (|24|) . In Sec. I11B1 we examine the 
validity of the next-to-leading order reduction formula using numerical calculation. In Sec. 
1111 Al we derive the reduction formula in the case of two resonances. The reduction formula 
Eq. (JilJ), which is the main result of this section, is expressed without the jump probability. 
In Sec. IIII Bl we numerically verify the validity of the reduction formula obtained in Sec. 
IIII Al comparing to that expressed by the jump probability. In Sec. IIV[ we summarize the 
results obtained in the present paper. 



II. REDUCTION FORMULA FOR ONE RESONANCE 

We derive the next-to-leading order reduction formula in the case of one resonance in 
this section. While the formula is suitable for investigation of the solar neutrino problem, it 
would be applicable for other cases. We then verify its validity using numerical calculation. 



A. Derivation of the next-to-leading order reduction formula 



We derive the next-to-leading order reduction formula in the presence of one resonance. 
Neutrino propagation in matter for three flavors is governed by 
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(3) 



in the base of the weak eigenstate <p a (a = e, n, r). Here, A(t) = \^2GpN e (t) is the matter 
induced mass of the electron neutrino with N e (t) being the electron number density. The 
quantities A i2 = Srn ^~^- m i) anc [ a 13 = ■ 5m i3(-ra 3 - m i) , are given in terms of the neutrino 
masses irti (i = 1, 2, 3) and the neutrino energy E. Throughout the present paper, we assume 
a hierarchical relation < Ax2 <C A13. The mixing matrix V is parametrized as 
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where Sy and represent sinfly and cos 0^, respectively. Since we discuss only the survival 
rate of the electron neutrino, which does not depend on #23 and 5, we set them to zero in 
the present paper. 

For our purpose, it is convenient to work in the base £ = Vi 3 V 23 (p, where the time evolution 
equation (jSJ) is rewritten as 
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Under the hierarchy A 12 , A(x) <C A 13 , the second matrix in the right-hand side (RHS) may 
be treated as a perturbation, while the 2x2 submatrix in the first matrix should not be 
dealt with as a small perturbation, because neglect of the submatrix causes a degeneracy of 
eigenvalues. Thus, neglecting the second matrix, we obtain a block diagonalized equation, 
and the leading order reduction formula Eq. (JTJ) is easily obtained. In this approximation, 
the discarded terms have the magnitude of order 0(Asi 3 Ci 3 ), which handles the error of this 
formula. We then derive a next-to-leading order reduction formula to reduce the magnitude 
of the error. 



For this purpose, let us rewrite Eq. (JJJ) as 
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In order to incorporate the effect of the perturbation v4si3Ci3, we then diagonalize the second 
matrix by moving to a new base 77 = VJ£, with V e being a time dependent unitary matrix, 

1 cose(t) sine(t) ^ 
K = 1 (9) 
y — sine(t) cose(t) j 



The diagonalization goes as 

1 A l2 s 2 12 + Ac 2 n As 13 c 13 \ I k_ N 
V} V e = 

^ ^si 3 ci 3 Ai 3 + As 2 13 ) \ k+ j 

where the two eigenvalues are given as 
1 



(10) 
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A + Ai 3 + s( 2 A 12 ± J (A + Ai 3 + S 12 A 12 ) 2 -4{Ai 3 A 12^12 + A(cf 3 Ai 3 + s 2 13 s 2 12 A 12 )} 



and the angle e satisfies the relation 

tan2e(t) 



A(t) sin 2^ 
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A 13 - sf 2 A 12 - A(t) cos 2^i3 ' 
Although the hierarchy A12, A <C A 13 implies that the time dependent angle is small, 

A sin 2#i3 
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the angle plays an important role in the improvement of the reduction formula. While the 
second matrix on the RHS of Eq. (jHJ) is diagonalized, we have additional terms in the time 
evolution equation in the base of 77. First, since the base 77 depends on time, the left-hand 
side of Eq. (jHj) yields the following extra contribution: 
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This term has a magnitude of the following order: 
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(15) 



where L denotes the typical length in which the matter effect A changes. In the case of the 
sun, the typical length corresponds to the scale height, L = R s ~ j^. 

Second, the first matrix in Eq. (JHJ is also slightly altered due to the change of base to 
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Ignoring the small extra off-diagonal terms, we approximately obtain the propagation 
equation in the time dependent base rj: 

( k_ A 12 s 12 c 12 + 0(A 12 sm29 12 e 2 ) O(^m) \ 
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where 



A e ff = k- - A l2 s 2 l2 . 



(19) 



This effective matter effect is well approximated under the hierarchy A\ 2 ,A(x) <C A13 as 

sin 2 26 13 A 2 



Aeff -> C 2 l3 A - 
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(20) 



and reduces to A e ff used in the leading formula in the limit — > 0. The off-diagonal terms 
discarded in Eq. |T8j) are smaller than the second matrix of Eq. ((ZJ) by the order Al2 ^° 3 2ei2 
in the case 

A sin 26>i3 



LA 



< A10 sin 26*1 2 e 



13 



which can be rewritten as 



E < LSm 2 l2 sin 29 



7 12- 



(21) 



(22) 



This condition is well satisfied for the solar neutrinos, since the parameters are E < 
10[MeV],L = R s ~ 3 x 10 14 [eV _1 ], and 5m 2 2 ~ l(T 4 [eV 2 ]. We also discard the element 
of 0(Ai 2 sin2#i 2 e 2 ) in the 2x2 submatrix because of its smallness. Thus, the order of the 
neglected matrix elements in Eq. (fTHJ) is much reduced, compared with that in the lead- 
ing order formula, and a considerable improvement is expected in the next-to-leading order 
formula obtained below. 
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The survival rate for the electron neutrino is calculated, by solving Eq. (jl8j) . to be 
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(23) 



Here, R 2x2 (t, 0) is the resolvent matrix for the two-flavor neutrino oscillation with the matter 
effect A e ff and mass squared difference A 12 . Thus we get the next-to- leading order reduction 
formula 



S 3v {v e 



cos 2 9 13 cos 2 (0 13 + e(0))S 2v>eff (A eff , A 12 ) + sin 2 9 13 sin 2 (0 13 + e(0)). (24) 



Here, S2 V)e ff(A e ff, A 12 ) is the survival rate for the electron neutrino in the effective two- 
flavor system with the matter effect A e ff and mass squared difference Ai 2 . This reduction 
formula has a twofold improvement compared with the leading order formula Eq. n(JIJ; the 
matter effect A e ff is replaced by A e ff and the angle is corrected from 9\ 3 to 9\ 3 + e at the 
production point of the neutrino, t = j^. 

This next-to-leading formula is valid no matter whether or not the resonance due to 
the matter effect occurs. We can, however, further simplify this relation if the effective 
survival rate S 2Ute ff(A e ff, Ai 2 ) is represented using the jump probability between the two 
mass eigenstates at the resonance point P c (A e ff, 9, A), which is almost the same as that for 
A e ff. The two-flavor survival rate is represented using the jump probability as 



S 2 u,eff = ^ + |^ - P c (A eff ,9i 2 , A i2 )| cos 29 m0 cos 26>i 2 

where 8 mQ is the mixing angle in the matter at the neutrino production point, 

A 12 cos20 12 -i e/ /(0) 
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We say that the resonance is complete in the case that the electron neutrino is produced far 
above the resonance point. In this case, the survival rate is represented by a jump probability 
as in Eq. ([25)1 . On the other hand, we say that the resonance is incomplete in the case that 



8 



the neutrino is produced too near the resonance point. In this case, we cannot represent the 
survival rate by the jump probability. The survival rate for an electron neutrino produced 
at a point with infinite matter effect, A e ff = c( 3 A — > oo, is represented by 

S 2 u,ef f = \-{\- Pc(A eff ,9 12 , A 12 )}cos20 12 . (27) 

For the case P c (A e ff, # 12 , A 12 ) — P c (A e ff, 9i 2 , A 12 ), the survival rate S 2u ,eff is rewritten as 

5W,e// = - ~ (s2u,eff — r) cos26? m0 - (28) 



2 V ' JJ 2 

Substituting this relation into Eq. (|2^j). we obtain the simplified formula 

S 3u {"e -> f e ) = cos 2 6» 13 cos 2 (6 l i3 + e) |i - (s 2 „, e ff - ^) cos20 mO } + sin 2 ^ sin 2 (0 13 + e). 

(29) 

In some cases, the survival rate S 2l/fi ff is known exactly, and the simplified formula can be 
represented in completely analytic form. For example, when the matter density distribution 
is of the exponential type, A = v4oexp(— ^), which we adopt in the numerical analysis in 
the next section, the survival rate is represented analytically using the jump probability, 

= exp(2 7 rAr cosg 2 )-l _ 
cK ' ' ; exp(27rAr ) - 1 V ; 

Using Eqs. (HH , flU} , flUJ) , (JZZj) , (EHJ , © , the survival rate for the three-flavor neu- 
trino oscillation is represented analytically under the condition P c (A e ff, 9± 2 , Ai 2 ) ~ 
P c (A e ff, $i 2 , A 12 ), in this case. We notice that this simplification is valid only in the case 
that the survival rate can be represented by the jump probability. It cannot be used in the 
case of incomplete resonance as is seen from the numerical analysis in the next subsection. 

B. Numerical confirmation of the validity of the next-to- leading order reduction 
formula 

In the present subsection, we numerically examine the precision of the next-to-leading 
order reduction formula. As an example, we use an exponential type electron density distri- 
bution N(r) = 245 N a exp(— 10.54^), which is a good approximation for the solar neutrino 



propagation in most regions 



•0 ' 



24] 39]. Here, Na is Avogadro's number. We consider the mass 



squared differences 1(T 10 < 5mj 2 < l(T 4 [eV 2 ], 1.5 x 1(T 3 < 5mj 3 < 4 x l(T 3 [eV 2 ] and the 



whole range of mixing angles < 612 < f , < #i 3 < |. We average over the survival rate 
with respect to the final time, and do not do it with respect to the initial time, since we do 
not use averaging in deriving Eq. (J2IJ). 

Since we have a vast parameter region, we do not investigate the whole parameter region. 
Instead, we investigate the parameter region where the error of the reduction formula is 
expected to be large, since we would like to be conservative concerning the precision of the 
reduction formula. We first fix the values of the mass squared differences as 5m\ 2 = 10 _4 [eV 2 ] 
and Sm1 3 = 1.5 x 10~ 3 [eV 2 ], where the hierarchy between the mass differences is minimum 
in the region we consider and the error is expected to be maximally enhanced. 

We show the error of the reduction formula, the deviation from the exact result Srf — 
Sexact, as a function of #i 3 in Fig. Ufa). The result derived from the leading formula Eq. 

is known to have the largest error of about 0.07. The next-to-leading order formula Eq. 
(j24j) without the replacement A e ff — > A e ff has an error of about 0.01. The next-to-leading 
order formula Eq. P^j) with the replacement A e ff A e ff has the smallest error, of order 
10~ 4 . We show the same figure for smaller value of #12 in Fig. ^b). The magnitudes of 
each error are almost the same as those of Fig. E^a). Theyarise, however, at relatively small 
value of #13 in this case. Generally, the errors are small for small sin 2#i3 as is expected from 
the fact that the angle correction e is proportional to sin 2#i3 (Eq. ([THjO and the correction 
to the matter effect is proportional to sin 2 2#i3 (Eq. (J2"U|) \ 

We show the results of the error of the reduction formula as a function of #12 for sin 2 2# 13 = 
1 in Fig. Efa), where the error is expected to be large. The leading order formula Eq. (P) has 
the largest error of about 0.06. The next-to- leading order reduction formula Eq. (|2"^j) without 
the replacement A e ff — > A e ff has an error around 0.01. The next-to-leading order reduction 
formula Eq. (|2*^|) with the replacement A e ff —>■ A e ff has an error of order 10~ 4 . Although 
the error of the leading order approximation vanishes for extremely large or small values of 
6\2i this is accidental. Since the survival rate S2 U , e ff is apparently unity for extremely small 
or large 612, the error comes from the angle correction e in this region. The correction of 
Sz„ due to the mixing angle correction is estimated to be 



sin 2 26 13 A(0) 



S 13 - C 13 S 2u,eff) (31) 



2A 13 

from Eqs. (JTHJ) , (J2H) • For #13 = | and S^.e// = 1, this quantity accidentally vanishes. This 
error therefore remains for other values of 613 as shown in Fig. Efb), which is the same figure 
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as Fig. |2fa) except for #i 3 = ~. 

We next confirm that the error tends to be smaller for a larger hierarchy, i.e., smaller 
5mf 2 and larger 5mf 3 , as is suggested by the fact that the errors are handled by the relative 
importance of the neglected off-diagonal elements compared with the dominant matrix ele- 
ment Ai 3 , i.e., A(0 )f 2ei3 for the leading order formula and Al2 T 20126 ~ Al2Sin2e f A 1 (0)sin2ei3 
for the next-to-leading order formula. We first show the 8m\ 2 dependence of the error in Fig. 
El The leading formula Eq. (0) has the largest error as in the above cases. There remains 
a finite error, even for extremely small values of 5m\ 2 . The errors for the next-to- leading 
order approximations tend to vanish for smaller values of 8m\ 2 . This is because the error 
is handled by Al2 sin2e ^4(°) hm2g i3 anc ^ SU gg es t s that the correction to the matter effect Eq. 
(|19j) is negligible in this case. Therefore, the correction to the matter effect will be negligible 
for the LOW and VO solutions of the solar neutrino oscillation. 

We next show the dependence of the error on the mass difference Smf 3 in Fig. 0] The 
leading formula Eq. ([TJ again has the largest error. The error tends to be reduced for larger 
values of 5m 2 3 , as is expected. 

All of the above numerical results strongly suggest that the reduction formula at the 
leading order potentially has an error around 0.1. This corresponds to the fact that the 
error due to the neglect of the mixing angle correction e(0) is estimated as 

sin 2 2fl 13 A(0) 2 2 - 

2^ Oia - <haS%>,eff) % 0.1, (32) 

from Eqs. (|13 )) . (j24j) . The numerical calculations also suggest that the next-to-leading order 
correction drastically improves the reduction formula and the error is reduced to be of order 
10~ 4 . These results indicate that the leading order reduction formula is pretty good for a 
rough estimation of the allowed parameter region and the next-to-leading order formula is 
necessary for its precise determination. A detailed analysis of the allowed parameter region 
from the solar neutrino data is not with in the scope of the present paper. 

Finally, we examine how precise the simplified reduction formula is j^j]. Since this 
formula is derived under the condition that the survival rate can be expressed by the jump 
probability, it will not be precise when the resonance is incomplete. For the case of the 
above example, N(r) = 245iVAexp(— 10.54-^-), the resonance is incomplete for neutrinos 
which have smaller energies ~ 1 [MeV]. The simplified reduction formula therefore is not 
expected to be valid in this case. On the other hand, the simplified formula is expected to be 
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FIG. 1: The error of the reduction formula is shown for the parameters written in the figure as a 
function of #13 in (a). The result using the leading formula Eq. ([T|). denoted by "LRF", has the 
largest error of about 0.07. The next-to-leading order formula Eq. 1)24|) without the replacement 
A e ff — * A e ff, denoted by "NLRF"', has an error of about 0.01. The next-to-leading order formula 
Eq. (|24|) with A e ffi denoted by "NLRF", has the smallest error of order 10 -4 . Larger errors arise 
for larger sin 26*13. (b) The same (a) for smaller value of #12. The magnitudes of each error are 
almost the same as those of (a). They arise, however, at relatively small values of #13. 
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FIG. 2: (a) The error of the reduction formula for the parameters written in the figure as a function 
of 9\2- The leading formula Eq. (JTJ) has the largest error of about 0.06. The next-to-leading order 
reduction formula Eq. (|24j) without the replacement A e ff — ► A e ff has error around 0.01. The 
next-to-leading order reduction formula Eq. Q24|) with the replacement A e ff — > A e ff has an error 
of only order 10~ 4 . Although the error of the leading order approximation vanishes for extremely 
large or small values of 612, this is accidental one. (b) Same as (a) for #13 = ~. Although the 
magnitudes of the errors are almost the same as in (a) there remains a finite error for extremely 
large or small values of 9±2. (b) shows the result for the same parameter set except for #13 = |. 
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FIG. 3: The error of the reduction formula is shown for the parameters written in the figure 
as a function of 5m\ 2 - (a) and (b) are the same except for the scale of the vertical axis. The 
leading formula Eq. Q has the largest error of about 0.05. The other lines, which have errors 
less than 0.01, are the next-to-leading order formula Eq. ()24|) with and without the replacement 



A, 



eff 



A e ff. For the leading formula, the error remains even for small <5mf 2 , although it tends 



to vanish for the next-to-leading order formula. 



14 



0.01 r 
a * 




-o.os 1 ' 1 ' ' 

1.5 2 2.5 3 3.5 4 

5m,y 10" 3 [eV 2 ] 

FIG. 4: The error of the reduction formula is shown for the parameters written in the figure as a 
function of 8m\%. The leading formula Eq. Q has the largest error of about 0.06. The other lines, 
which have errors less than 0.01, are the next-to-leading order formula Eq. (|24|). with and without 
the replacement A e tf — > A e ff. All errors tend to be reduced for larger <5mf 3 . 

precise for neutrinos that have larger energies ~ 10 [MeV], because in this case the resonance 
is expected to occur almost completely. We show in Fig. El the error of the reduction formula 
as a function of the neutrino energy E. One can observe that the simplified formula is not 
valid for smaller energies of the neutrino from this figure. The errors both of the leading 
and the next-to-leading order reduction formulas are also shown in the figure. According to 
the above result, we learn that we should not use the simplified formula without ensuring 
that the resonance is complete, although the formula is attractive because of its simplicity. 

So far we have not restricted the values of the mixing angles, since our aim was to 
conservatively clarify the precision of the next-to-leading order reduction formula. However, 
there exist meaningful upper bounds on the mixing angle #13 from the reactor experiments 
CHOOZ j^jj and Palo Verde j^. The smaller the mixing angle #13 becomes, the better 
the precision of reduction formulas is expected to be. Thus we have calculated the errors of 
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FIG. 5: The error of the reduction formula is shown for the parameters written in the figure as 
a function of the neutrino energy E. The error of the simplified formula is very large for smaller 
values of the neutrino energies. On the other hand, the leading order reduction formula has an 
error around 0.06. The next-to-leading order reduction formula has an error of order 10~ 4 . 

the reduction formulas for sin 2 2# 13 = 0.12, the upper bound corresponding to the best fit 
value 5m1 3 = 3 x 10~ 3 [eV 2 ] to account for the atmospheric neutrino data |33j]. The results 
are shown in Fig. H (a), as functions of the remaining mixing angle tan 2 #i2. Also shown 
by a shaded area is the region allowed at 90% C.L. by the LMA MSW solution Q- We 
learn from this figure that the leading order formula has an error up to 1% or so for some 
values of tan 2 #i2, while the error of the next-to-leading order formula Eq.(|24|) is essentially 
negligible. The error of the leading-order formula is, however, less than 0.2% if we remain 
in the shaded region. In Fig. EJb), we have also shown the errors of the reduction formulas 
for sin 2 2#i3 = 0.28, the upper bound corresponding to the value 5m\ 



l 13 



1.5 x 10~ 3 [eV 2 l, 

n 

the lowest mass squared difference to account for the atmospheric neutrino data |3J3] at 90% 



C.L. We now learn that the error of the leading order formula is enhanced by both smaller 
<5m 2 3 and larger # 13 ; the error reaches 5% or so. Even for the shaded region, the error can 
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be up to 1%. 

From these analyses, we can say that for the values of #13 implied by the reactor exper- 
iments the error of the leading order formula can be rather small, while that of the 
next-to-leading order formula is completely negligible. If we further impose the condition 
suggested by the LMA MSW solution, the error is even smaller, i.e., at most 1%. Such an 
error, however, will be problematic for precision tests of neutrino experiments, such as the 
(indirect) search for CP violation 0,0, HQ, which require the precise determination 



of the size of the unitarity triangle with an accuracy of a few percent [23], and therefore 
a better precision of the reduction formula itself. We hope that the next-to-leading order 
reduction formula proposed here, being a simple formula to use, will be useful for study of 
the precision tests of neutrino experiments. 



III. REDUCTION FORMULA IN THE CASE OF TWO RESONANCES 



We derive a reduction formula when v e experiences two successive resonances in Sec. 
IIII Al The expression we derive is slightly different from previously proposed ones and is 
applicable for wider situations. This formula is relevant for investigation into the supernova 
neutrino data or the h ypo thetial very high energy solar neutrino data due to the annihilation 



of WIMPs 125 
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We then verify the validity of the reduction formula using 



the numerical calculation in Sec IIII Bl 



A. Derivation of the reduction formula in the case of two resonances 

Suppose v e is produced at time and detected at time t, going through "higher" and 
"lower" resonances, caused by matching the matter effect with A13 and A12, respectively. 
We divide the time interval into [0,^] and \pM,t] (0 < tu < t), where the conditions 
A 12 <C A 13 , A and A 12 , A <C A 13 are met, i.e., the higher and lower resonances are operative. 
The intermediate time tu is chosen so that A i2 <C A(tu) <C A13. 

We first consider the time range where the higher resonance occurs, i.e., A12 <C A13, A(t). 
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FIG. 6: The error of the reduction formula for the realistic parameters (a) &m\ 3 = 3.0 x 
l(r 3 [eV 2 ],sin 2 2f9i3 = 0.12 and (b) 5m 2 l3 = 1.5 x 10 _3 [eV 2 ], sin 2 20 13 = 0.28. The error of the lead- 
ing formula Eq. Q is shown as " LRF" . The errors of the next-to-leading order formula Eq. (|24jl 
without and with the replacement A e ff — > A e ff are shown as "NLRF"' and "NRF", respectively. 
The allowed region from the solar neutrino data is shaded. The error for the case (b) is found to 
be larger than that of (a), and the errors are rather small for the shaded region in both cases. 
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uj = U(f), 



u 



A convenient base to describe this region is 

f -I q\ 
1 

V° l J 

where the time evolution equation Eq. (JSj) can be cast into 



v2 



23- 



(33) 



dt 





f° 1 






> 


< uv 

















{ oj 





(34) 



-S12C12C13A12 



S12C12S13A 



12 



-Si2C 12 Ci3A 12 S^ 2 C^ 3 A 12 + S^ 3 A 13 + Ait) S 13 Ci3A 13 — S^S^C^A]^ 

V S12C12S13A12 S13C13A13 - ^ 2 si3Ci 3 Ai2 s\ 2 s\ 3 A l2 + c 2 13 A 13 



uj (35) 



cos2# 12 A 12 











V 



s? 3 (A 13 - S 2 2 A 12 ) + A(t) s 13 c 13 (A 13 - s? 2 A 12 ) 

c? 3 ( A i3 - s 2 12 A 12 ) ) 





Si 3 Ci 3 (Ai 3 - S 12 L^i2 

+s\ 2 A 12 Iuj 

<os2# 12 A 12 T 
M H (t) 

Here, we neglect the term proportional to the unit matrix /. 
The resolvent matrix in this base R H is 



UJ 



UJ. 



(36) 
(37) 



R H {t,0) 



-i cos 20i2 Ai2t 



V 







(38) 



T 

>(t,o) 

where R 2 ¥ x2 is the resolvent matrix for the effective two-flavor neutrino system with matter 
effect A(t), mass squared difference A13 — s^ 2 A 12 , and mixing angle 9 13 . 

We next consider the time range, [£m,£], where the lower resonance occurs, i.e., 
A12, A(t) <C A13. A convenient base in this region is the time dependent base n(t) = Z(t)<f), 
by which the Hamiltonian for fixed t is diagonalized, 



(&4 < k 2 < k 3 







(0 ] 




' Ait) > 


- 









\ 


Z(t)< 


V 


A12 









> z\t) = 





h(t) 









{ A 13 J 










I 





h(t) ) 



19 



Because of the hierarchy A 12, A <C A 13, the Hamiltonian and therefore the resolvent are 
approximately block diagonalized: 



R 2x2 (t, tu) 
T 





e -if tM Mt>)dt> 



(40) 



Since the lower resonance is nearly complete, i.e., A 12 <C A{t M ) <C A 13 , the 2x2 resolvent 
matrix can be generally written as 



1 



(41) 



VI - P L e ia y/PleW 
yfFle^ -VI - P L e~ i{ - a - p -^ 

where ti is some time far after the lower resonance, i.e., A(ti) <C A 12 , and is the jump 
probability between adiabatic states with respect to the lower resonance. 
Using these resolvent matrices, the survival rate of v e is given as 

\ 



S 



1 ouW RUMa/) ° 

^ T eM-HLHtW)) 



Z(t 



-i cos 2612 Ai2*m 



0^ 



M) 



u 



^-2x2(^,0) 







(42) 



Here, the matrix Z{tu)U^ is approximately written as 

( I N 

C13 -S13 



^JVf)^ 



1 0^ 



V 3m (t M ) 



(43) 



\0 S13 C13 J 

where V 3m (t) is the matrix that diagonalizes the 2x2 mass matrix M# as 

/ 



Vi m (t)M H (t)V 3m {t) 



h(t) 
V h(t) 

From Eq. g^) and Eq. and [7(1, 0, 0) T = (0, 1, 0) T , 

/ 

/ RL-ft.**/) \ 

S ~ (10 ^(*) 



(44) 



R 2x2 (f, £m) 

T exp(-ij* fc(f)dO 



V 



V 3m (tM)Rjx2(£M, 0) 



v ^ 



(45) 
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In order to express the survival rate in terms of the two-flavor survival rate, we set 



V 



($M ) R-2X 2 (*M , 0) 



(46) 



The survival rate for the two-flavor system with respect to the higher resonance is expressed 
by these quantities as 



>2x2 



1 )R& 2 (t,0) 



i o ) v 3m (t)vL(t)Rf X 2(*,o) 



V 



1 V 3m (t) 



/ -i r* k 2 (t')df 



V 



o 



-i f* k 3 (t')df 

e Jt M 



a\ 2 cos 2 6 13 + \b\ A sm A 6 13 . 



In Eq. (J5UJ) . we averaged with respect to the final time t and used 



V 3m (i) ~ V 3m (t 



Cl3 

\ — Si3 C13 

Using the unitarity condition \a\ 2 + \b\ 2 = 1, we get 

l fl |2 _ ^2x2 ~ sm2 ^13 |^|2 _ cos2 #13 - »S'2x2 



cos 20 13 

The survival rate is thus expressed as 



cos 20 



13 



5 



(l oW(*) 



/ L / 



a 
6 



a| 2 cos 2 6 13 {P L cos 2 12 + (1 - P L ) sin 2 12 } + |6| 2 sin 2 13 . 



(47) 
(48) 

(49) 
(50) 



(51) 



(52) 
(53) 



In Eq. we averaged the survival rate with respect to the final time and used Z\t) = V. 
Using Eq. (|5Tj). we get the reduction formula, which is our final result, 



^2x2 — s i n2 ^13 2/1 , cos2 ^13 ~ 5 , 2x2 • 2 



cos 20 13 



sin 0i 3 . 



(54) 
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Here, 5 , 2x2,m/( = Plcos 2 9\ 2 + (1 — Pl) sin 2 #i 2 ) is the two-flavor survival rate for the lower 
resonance in the case that the electron neutrino is produced at a point with A —>■ 00. We 
verify the validity of this formula by a numerical method in the next subsection. 
This reduction formula coincides with the reduction formula obtained in Ref. 

S = P 2 H cos 2 9 13 S^ mf + (1 - Pi) sin 2 9 13 , (55) 
(P 2 H = ^ + ^-P H )cos29 13m y 

provided the higher resonance is complete and the two-flavor survival rate $1x2 can be 
written in terms of the jump probability Ph as 

S&2 = \ + {\~Ph) cos2# 13m cos 20 13 . (56) 

This reduction formula Eq. (}55|) is the same as the reduction formula Eq. (J2j) obtained 
in Ref. p, if # 13 is set to be | and S 2X 2,mf * s further expressed by the jump probability 
Pl as S , 2x2,in/ = \ ~ (| — -Pl) c °s26 l i2. Although the reduction formulas Eq. (JH3|) and 
Eq. are the same when the higher resonance is complete, they are different when the 
higher resonance is incomplete [4^. We also confirm the difference numerically in the next 
subsection. 



B. Numerical confirmation of the validity of the reduction formula for the case of 
two resonances 

In the present subsection, we examine the precision of the reduction formula for the case 
of two resonances Eq. (JH3J) using a numerical calculation. As a typical example, we use the 
same electron density distribution N(r) = 245iV4 exp(— 10.54-^-) as in Sec. Ill Bl Since we 
consider the case where two resonances occur, we assume that the energy of the produced 
electron neutrino is very high compared to that of the solar neutrino, E > 100 [MeV]. This 
situation is quite similar to that considered in Ref. To get conservative results, we 

take the mass squared differences as 5m\ 2 = 10 _4 [eV 2 ] and <5m 2 3 = 1.5 x 10~ 3 [eV 2 ], where 
the hierarchy between the mass squared differences is mildest in the experimentally allowed 
region and the error is expected to be enhanced maximally. 

We show the error of the reduction formula Eq. (|54|) in Fig. [7|for various values of the 
electron energy E = 0.1, 1, 10, 100 [GeV]. We observe the largest error for the lowest energy 
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E = 0.1 [GeV]. Even in this case, however, the error is less than 0.005. For higher energy 
neutrinos, the reduction formula is more accurate. This reduction formula is apparently 
accurate enough for the investigation into supernova neutrinos and the hypothetical O(GeV) 
solar neutrinos 25, 26, 27, 28,Q,Q- 

We finally compare the reduction formula obtained here [Eq. (J54j) ] with that obtained in 
Ref. |30|. We compare the errors as a function of the initial electron density iV(0) = 245GA0i 
in Fig. |HJ The reduction formula Eq. ()54|) is extremely accurate for all values of 0. While 
the reduction formula Eq. ()55|) is better than Eq. (j2J), it is not accurate for smaller values 
of G where the higher resonance becomes incomplete. The error, however, is small enough 
in the parameter range considered in Ref. jscj for their purpose. This result shows that the 
reduction formula Eq. (|54|) can be safely used even in the region where the higher resonance 
is incomplete. 

IV. SUMMARY 



e 



In the present paper, we derive the next-to-leading order reduction formula for the v { 
survival rate Eq. (J2"4~)) from the three-flavor neutrino oscillation to the two-flavor one in 
the case when there is only one resonance, as in the ordinary solar neutrino oscillation. 
Together with an analytic argument, we numerically verify the accuracy of the reduction 
formula, leaving the mixing angles free for generality. While we find that the leading order 
reduction formula Eq. (Q) is accurate enough for a rough estimation, the next-to-leading 
order reduction formula is extremely accurate and adequate for precision tests of neutrino 
oscillations. Next, we study the accuracy of the reduction formulas in a realistic case, i.e., 
taking into account the current upper bound on # 13 . We find that the largest error of the 
leading order reduction formula is about 1% or so while the error of the next-to-leading order 
reduction formula is negligible. We thus point out that this precise next-to-leading order 



formula will be useful 
study of CP violation 



or precision tests of neutrino oscillations, for example the (indirect) 
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23j . We also verify the accuracy of the reduction formula 
written using the jump probability. This formula is accurate when the resonance is complete, 
i.e., for high energy neutrinos, although it is not valid when the resonance is incomplete, 
i.e., for low energy neutrinos. 

We also derive the reduction formula Eq. ()54j) in the case of two resonances as in the 
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FIG. 7: The error of the reduction formula Eq. (|54[) is shown for the parameters written in the 
figure as a function of 812 and #13 for various values of the neutrino energy E = 0.1,1,10,100 
[GeV]. The largest error occurs in the case of the lowest energy E = 0.1 [GeV]. Even in this case, 
the error has magnitude less than 0.005. 
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FIG. 8: The errors of the reduction formulas Eq. (|54jl. Eg. (j55j) . and Eg . ((2j) as a function of the ini- 
tial electron density. The horizontal axis is which is defined as N(r) = 2450 Na exp(— 10.54-^). 
The lines (1), (2) and (3) correspond to Egs. (|54j). ©, and (|55|l. respectively. The reduction 
formula Eg. (|54|1 is extremely accurate for all values of 0. While the reduction formula Eg. (|55|) 
is better than Eg. (J2J), neither is accurate for smaller values of 0. 

oscillations of the supernova neutrinos and the hypothetical high energy O(GeV) solar neu- 
trinos due to the annihilation of WIMPs. We numerically verify that it is quite accurate and 
applicable for any parameter region. We then compare it to the reduction formulas obtained 
in Refs. Although the previously obtained formulas are valid only in case the higher 

resonance is complete, the formula obtained here is valid not only for the complete case but 
also for the incomplete case. 
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